Pressure-dependence of Curie temperature and resistivity in complex Heusler alloys 
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Using first-principles electronic structure calculations, we have studied the dependence of the 
Curie temperature on external hydrostatic pressure for random Ni2MnSn Heusler alloys doped with 
Cu- and Pd-atoms, over the entire range of dopant concentrations. The Curie temperatures are 
calculated by applying random-phase approximation to the Heisenberg Hamiltonian whose param- 
eters are determined using the linear response and multiple scattering methods, based on density- 
functional theory. In (Nii-^Pd^MnSn alloys the Curie temperature is found to increase with 
applied pressure over the whole concentration range. The crossover from the increase to the de- 
crease of the Curie temperature with pressure takes place for Cu-concentrations larger than about 
70% in (Nii_3;,Cua;)2MnSn Heusler alloys. The results for the reference Ni2MnSn Heusler alloy agree 
well with a previous theoretical study of Sasioglu et al. and also reasonably well with available ex- 
perimental data. Results for the spin-disorder induced part of the resistivity in (Nii-^Pd^MnSn 
Heusler alloys, calculated by using the disordered local moment model, are also presented. Fi- 
nally, a qualitative understanding of the results, based on Anderson's superexchange interaction 
and Steam's model of the indirect exchange interaction between localized and itinerant d-electrons, 
is provided. 

PACS numbers: 71.55.Ak,72.25.Ba,75.10.Hk,75.30.Et 



I. INTRODUCTION 



High pressure studies form an important area of re- 
search in solid state physics. Over the years such studies 
have provided useful insight into our understanding of the 
physical properties of solids*^ Pressure has basically two 
effects on the electrons in solids: increased kinetic energy 
and accompanying change in the effects of Coulomb in- 
teraction among the electrons. In typical band structure 
calculations the former effect is captured via increased 
overlap of the basis orbitals, leading to increased hopping 
matrix elements and band broadening. The changes in 
the correlation effects are captured via self-consistency 
of charge and potential, as dictated by the density func- 
tional theory and its variants (e.g. generalized gradient 
approximation (GGA), local density approximation with 
on-site Coulomb interaction (LDA+U)). Both of these 
affect magnetic properties of materials such as exchange 
interaction, local moments, Curie temperature and resis- 
tivity due to magnetic scattering. Of particular impor- 
tance for metallic magnets are the change in the Fermi 
surface and the hybridization between different orbitals, 
which influences the itinerant vs. localized nature of the 
charge carriers. It is a useful exercise to explore how the 
existing theories of electronic and magnetic structure cal- 
culation fare in describing the pressure variation of these 
physical properties. In the present paper we study the 
pressure-dependence of the Curie temperature and the 



resistivity of some quaternary Heusler alloys, which form 
an important class of magnetic materials, with poten- 
tial industrial/technological applications^— The present 
work is a natural continuation of our previous paper— in 
which we presented an extensive study of the magnetic 
and transport properties of these alloys at ambient pres- 
sure. 

There exists a number of experimental studies of the 
variation of the Curie temperature (T c ) under pressure 
for elemental ferromagnets, transition-metal alloys^ 
as well as Heusler alloys^— including Co2TiAljiS 
Ni2MnSn^iii Pd2MnSn and some others^ Recently 
diluted magnetic semiconductors^ and some random 
Heusler alloys^ have also been studied. However, 
the corresponding experimental studies of the pressure- 
dependence of the resistivity are rare: one example of 
such a study is the work by Austin and Mishra^ 

Several studies of the ambient/equilibrium electronic 
and magnetic properties of such alloys have appeared 
recently Theoretical studies of the pressure- 

dependence of the Curie temperature of elemental transi- 
tion metal ferromagnets include model-based studie s 19 ' 20 
as well as recent first-principles calculations! 21 ' 22 Turck 
et al£2- have studied the pressure-dependence of the Curie 
temperature in hep Gd and found a strong dependence 
of T c on the c/a ratio, which, they suggested, can be 
related to the measured T c of thick epitaxial Gd(0001) 
films on various transition-metal substrates. Similar cal- 
culations have also been carried out for the intcrmctallic 
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compound GdAl^ and for Cr-based compounds in the 
zinc blende structured A systematic theoretical study 
of the pressure-dependence of T c in Heuslcr alloys was 
carried out recently by Sasioglu et al.^ where the ex- 
perimentally observed increase of T c with pressure^^ in 
Ni 2 MnSn was reproduced correctly, albeit only qualita- 
tively so. There exists some other studies of Hcusler al- 
loys in which the dependence of the exchange integrals 
on volume contraction/expansio n 17 i 27 ~ — was explored. 

The particular Heuslcr systems considered in this work 
are the (Nii-^Cu^MnSn and (Nii_ x ,Pd x ) 2 MnSn al- 
loys. We employ the coherent potential approximation 
(CPA) to treat the disorder in the Ni-sublattice due to 
doping with Cu- or Pd-atoms. It is well known that CPA 
neglects the effects of short-range order, which are often 
present in random alloys. However, for the cases under 
study, such effects should be weak, as we need to describe 
the magnetic interactions on the non-random Mn sublat- 
tice, which is only indirectly influenced by the random- 
ness of the transition-metal sublattice via hybridization. 
The most important effect of disorder is, therefore, the 
change in carrier concentration, which is well-described 
by the CPA. In a recent paper^S it was demonstrated 
that the CPA provides a reasonably good description of 
the electronic structure for the case of disorder on the 
magnetic sublattice as well. 



II. FORMALISM 

The electronic structure calculations are performed 
employing the tight-binding linear muffin-tin orbital 
(TB-LMTO) basis^i and the density functional theory 
(DFT). The local spin-density approximation (LSDA) for 
the exchange-correlation part of the potential is used. 
The effect of substitutional disorder among (Ni,Cu) or 
(Ni,Pd) atoms is described by the coherent potential ap- 
proximation (CPA) as formulated in the framework of 
the TB-LMTO Green function method^ The calcula- 
tions employ an s,p, d, /-basis, the same atomic radii 
are adopted for all atoms, and the Vosko-Wilk-Nusair 
exchange-correlation potential 3 - 3 - is used. 

The exchange interactions and Curie temperatures are 
studied here by employing a classical Heisenberg Hamil- 
tonian 

H = -^2j ij e i -e j , (1) 

i¥=3 

where denotes the exchange integral between Mn 
atoms at sites i and j, and and are unit vectors 
in the directions of the local magnetization on sites i and 
j, respectively. The parameters of the Heisenberg Hamil- 
tonian are obtained using a two-step approach] 34 ! 35 The 
reference state for our calculations is chosen to be the 
disordered local moment (DLM) stated which was suc- 
cessfully used in previous studies The exchange in- 
teractions Jij can be expressed in the framework of the 



TB-LMTO method as 

Jij = ^-lm f tr L {A 2 (z)4.(z)A,(z)^(z)}dz, (2) 

where Aj(z) characterizes the exchange splitting of the 
Mn atom at the site i, gfj (z) is the configurationally aver- 
aged Green function describing the motion of an electron 
with spin er, a =tj4-) between Mn sites i and j, and the 
integration is done over the contour C in the complex 
energy plane z which starts below the valence band and 
ends at the Fermi energy. Symbol tr^ denotes the trace 
over the basis orbitals of angular momentum symmetry 
L = (£,m). The effective exchange splitting A, (z) is de- 
fined in the TB-LMTO method in terms of the potential 
functions Pf(z) of Mn atoms as A !; (z) = Pj(z) - P^iz). 
Further details on the exchange interactions evaluated in 
the DLM reference state can be found in Refl37l 

Wc have found that the mapping can be improved (pro- 
viding a better agreement of calculated T c with the ex- 
periment) by including the electron correlation effect on 
Ni-sites in the framework of the LSDA+U method (see 
Fig. 1 in Ref|(|). The electron correlations on Pd-sites 
are found to have negligible influence on the calculated 
T c . We determine the Curie temperature corresponding 
to the Heisenberg Hamiltonian using the RPA 3 - 5 . and by 
including all exchange interactions up to ~ 4 lattice con- 
stants. The pressure is simulated by reduction of the 
lattice constants from their ambient pressure values, as 
reported in experiments (see Table 4 of Refl38l and Table 
2 of Refl39l) . In ReffTTI an empirical pressure- volume re- 
lation was given for Ni 2 MnSn. Although such a relation 
can be obtained entirely from theoretical calculations, in 
the present study, as in Refill, we simply reduce the 
lattice constant and estimate the corresponding pressure 
from the empirical pressure-volume relation. According 
to both the empirical pressure- volume relation 1 ^ and the- 
oretical calculations r& a 3 % reduction of the lattice con- 
stant corresponds roughly to a pressure of 16 GPa in 
Ni 2 MnSn. 

In Heuslcr alloys (Nii_ x ,T x ) 2 MnSn (T=Cu, Pd) the 
effect of disorder among the Ni- and T-atoms at the 
Fermi energy Ep is weak because the Ni- and T-states lie 
well below Ep^ Consequently, the corresponding resid- 
ual resistivity is small. The resistivity due to phoiion 
scattering 4 ^ is known to be small as well. The dom- 
inant contribution to the resistivity, thus, is from the 
spin-disorder scattering, which increases with temper- 
ature, reaching its maximum value at T c . In Ref|U 
we employed a simplified approach to estimate the 
temperature-dependent resistivity in (Nii_ x ,Pd x )2MnSn 
Heusler alloys. It was shown that the spin-disorder at 
and above T c can be described satisfactorily using the 
DLM model, and the corresponding resistivity can be 
described by the Kubo-Greenwood approach adapted to 
the TB-LMTO method.— We have used the same ap- 
proach to study the resistivity under pressure as well. 
For further details concerning computational techniques 
we refer the reader to our recent paper £ 
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The method we have used is based on mapping the 
(T = 0) total energy to an effective Heisenberg Hamil- 
tonian. In fact it is just the band energy, not the total 
energy, which is considered in the mapping, by appealing 
to the so-called 'magnetic force theorem' (see Ref. (|25f ) and 
references therein). For finite temperature (T ^ 0) cal- 
culations, as would be appropriate for the DLM reference 
state, one needs to consider the mapping of the free en- 
ergy F = U — TS, where U is the internal energy and S is 
the entropy. U should include the electronic and the av- 
erage vibrational energy and S should include both elec- 
tronic and vibrational entropy. Electronic energy should 
be calculated via finite temperature density functional 
theory. This free energy cannot, in principle, be mapped 
to an effective Heisenberg form, which includes only the 
electronic parameter (spin or magnetic moment). Such 
a scheme could be implemented for a 'supercell' calcu- 
lation. However, the advantage of being able to calcu- 
late exchange interactions up to a large distance with 
relative ease and reasonable accuracy and therefore pre- 
dict results and trends for a group of materials will be 
gone. For finite temperatures, both electronic excitations 
and lattice vibrations will influence magnetic properties. 
However, in the absence of an accurate scheme, we will 
refrain from speculating on these effects on our results. 



III. RESULTS AND DISCUSSION 

In this section we present results for the pressure- 
dependence of the Curie temperature and spin-disorder 
resistivity of (Nii-^T^MnSn (T=Cu, Pd) Heusler al- 
loys over a broad range of concentrations. 



A. Curie temperature under pressure: previous 
studies 

We start with a brief review of previous theoretical 
studies, both model-based and first-principles. Theoret- 
ical study of the pressure-dependence of T c on a model 
level has usually been based on a discussion of the indi- 
rect exchange interaction among spins^ or on the Hub- 
bard model of magnetising It is clear that the pressure- 
dependence of T c carries information about the volume- 
dependence of the electron wave functions^ (Bloch func- 
tions in the case of crystalline solids). In an RKKY- 
type indirect exchange interaction model, the volume- 
dependence of T c can be related to the volume depen- 
dencies of the density of states at the Fermi level N(Ep) 
and the exchange integral or the matrix element of the 
exchange interaction between the wave functions at the 
Fermi energyi^(henceforth referred to as the bare ex- 
change interaction J barc , see Eq.Q). The density of 
states at the Fermi level usually decreases with pres- 
sure (unless, under pressure, new bands start crossing the 
Fermi level), as the bands broaden due to increased over- 
lap between the orbitals centered on neighboring atoms. 



The bare exchange integral usually increases under pres- 
sure, due to increased overlap of the wave functions. The 
above two effects compete with each other in determining 
the final change of T c under pressure. 

In a description based on the Hubbard model, the band 
terms and the term involving the intra-atomic Coulomb 
interaction U are written separately. In this case it is pos- 
sible to divide the volume/pressure-dependence of T c into 
a term originating from the band broadening under pres- 
sure and another that relates to the volume-dependence 
of U . The importance of this latter term for the pressure- 
dependence of T c in Ni and Ni-Cu alloys was first pointed 
out by Lang and Ehrenreichi^ They argued that it is this 
term involving the intra-atomic Coulomb interaction that 
is responsible for the increase in T c under pressure (P) 
in Ni and the decrease of dT c /dP with the increase of 
Cu-concentration in Ni-Cu alloys. 

In ab initio calculations based on DFT, such as ours, 
the above effects are merged together and it is hard to 
delineate one from the other. The advantage, of course, is 
that the exchange interactions are realistically calculated 
for the system at hand, and no assumption related to the 
volume-dependence of any parameter needs to be made. 
Prototypical studies in this direction are Refs l2lll22l ad- 
dressing the pressure-dependence of the Curie temper- 
ature of bec iron. In particular, a very recent first- 
principles study-22 seems to agree with experiment, which 
gives essentially no pressure-dependence of T c . The au- 
thors describe this result as a balance between the band 
structure effect which reduces the magnetization and an 
increase of the bare exchange integrals with decreasing 
volume, as mentioned above. A somewhat different sit- 
uation is encountered for fcc-Ni. This metal is a strong 
ferromagnet, its magnetic moment depends on pressure 
only weakly and an increase of exchange integrals un- 
der pressure thus dominates. As a result, T c of fcc-Ni 
increases with pressure. In general both effects are rele- 
vant and the observed change of T c under pressure is a 
result of the competition between them. 

The above-mentioned competition of the two effects is 
illustrated in Fig. [TJ Here we plot the effective exchange 
interactions Jij ( Eq(2]) , which control directly the Curie 
temperature of the system (Fig. QJi), and the bare ex- 
change interactions Jy are , defined as 

= Jij/{MiMj), (3) 

as a function of the distance between the Mn atoms in 
Ni2MnSn. We consider the ambient pressure case, as 
well as elevated pressure simulated by reduction of the 
lattice constant. In Eq.© Mi denotes the size of the 
local magnetic moment of Mn atom at site i (Fig. [TJd). 
This definition is motivated by the presence of the ex- 
change splittings Ai(z) in Eq. @ which are roughly pro- 
portional to the moment magnitudes Mj. The increase 
of the bare exchange interactions with pressure is clearly 
seen for the first five neighbors (see Fig. [TJd) . It is also 
obvious that magnetic moments decrease with the pres- 
sure, as expected. These two effects, i.e., enhancement of 
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FIG. 1: The Mn-Mn exchange interactions in the Ni2MnSn 
Heusler alloy as functions of the relative interatomic distance 
d/a where a is the fee lattice parameter: (a) the effective inter- 
actions Jij calculated from Eq. ((2j), (b) the bare interactions 
j]j are defined by Eq. ([3| ■ The cases of ambient pressure (open 
circles) and the pressure corresponding to the 3% reduction 
of the lattice constant (filled circles) are shown. 



jbare anc j SU pp r ession of magnetic moments, both due to 
the volume decrease, compete with each other, resulting 
in smaller pressure- induced changes of the effective inter- 
actions Jij (Fig. [TJt) . The final result in the present case 
is an enhancement of T c under hydrostatic pressure. The 
same qualitative explanation was put forward by Cson- 
tos et al^- for the increase of T c with pressure in some 
diluted magnetic semiconductors. 



B. Calculated results 

In this section we present and compare results for 
(Ni, Pd) 2 MnSn and (Ni, Cu) 2 MnSn Heusler alloys, with 
Ni2MnSn considered as the reference case. 



1. Ni 2 MnSn and Pd 2 MnSn Heusler 



The Heusler alloy N^MnSn was studied theoretically 
in great detail in a recent paper ^ while corresponding 
experimental results can be found in 

RefsiandOl]. The 

experiments give the rate of increase of T c with pres- 
sure around 0.62 K/kbar— (for pressures up to 10 kbar 
or 1 GPa) and 0.744 K/kbar for higher pressures up to 
9 GPaAi The authors of the theoretical study2£ employ 
the frozen-magnon approach and inverse lattice Fourier 
transform to construct the effective Hciscnbcrg Hamil- 
tonian, which is then analyzed in the framework of the 
multisublattice MFA. They obtain an increase of T c by 
38 K for pressures around 16 GPa, an increase from the 
ambient value T c =362 K to T c =400 K. Our estimated val- 
ues are T c =322 (334) K for ambient pressure and T c =375 
(401) K (an increase of 53 (67) K for 3% reduction of the 
lattice constant, roughly corresponding to the same pres- 
sure ~ 16 GPa). The values in brackets correspond to 
the model assuming electron correlations on Ni-sites (the 
effective Hubbard J/ni chosen is 2 eV). Assuming linear 
relation between T c and pressure, the experimental in- 
crease is expected to be around 100—120 K. The theory 
thus correctly predicts a pressure- induced increase of T c , 
while the absolute increase is at least three times (Refl26l) 
to 1.5— 1.8-times (present study) too small. The exper- 
imental data for Pd2MnSn Heusler alloys that we have 
come across in the literature^ also indicate an increase of 
T c with pressure at a rate of 0.75 K/kbar, the sameii or 
almost the same^ as that for Ni2MnSn. We have correctly 
obtained an increase of T c for 3% reduction of the ambi- 
ent pressure lattice constant, as in the case of Ni2MnSn. 
Specifically, T c is increased by 40 K, a value smaller than 
but still comparable to that for Ni2MnSn and in reason- 
able agreement with experiment. We find that any addi- 
tional (i.e., in addition to that given by LSDA) Coulomb 
term U on Pd-sites has a negligible effect on T c or its 
pressure-dependence. 

The hockey-stick-like variation of the composition de- 
pendence of T c for Nii-zjPd^MnSn alloys at ambient 
pressure is captured in the RPA, but not in the MFA. The 
RPA is equivalent to including the reaction field (or the 
Onsager cavity field) effect in the calculation of T c (see, 
for example, Cyrol4^ and references therein). Our results 
point to the importance of including the reaction field ef- 
fect in reproducing the correct (experimentally observed) 
composition-dependence of T c in this alloy system. 

The temperature can influence the above results, as 
discussed at the end of section HH Increasing temper- 
ature leads to an expansion of the lattice, compensat- 
ing somewhat for the pressure-caused contraction. In 
essence, thus, the above reasoning should be applied to a 
net effective reduction of volume, considering both tem- 
perature and pressure effects. Lattice expansion effects 
and the effect of including the Fermi-Dirac distribution 
function in performing the exchange integrals in Eq. (2) 
have been discussed recently by Ailing et a/..— However, 
a complete T ^ calculation, as described in section [TTl 




Pd concentration (x) Cu concentration (x) 



FIG. 2: The concentration-dependence of Curie temperatures 
(RPA) for (Nii-^Pd^MnSn alloys for ambient pressure and 
for the pressure corresponding to the reduction of the alloy 
lattice constant by 3%. Model assumes electron correlations 
on Ni-sites treated in the framework of the LSDA+U method. 



FIG. 3: (Color online) The concentration-dependence of Curie 
temperatures (RPA) for (Nii-zjCu^MnSn alloys for ambi- 
ent pressure and for set of pressures corresponding to indi- 
cated reductions of the alloy lattice constant. Model assumes 
electron correlations on Ni-sites treated in the framework of 
the LSDA+U method. 



has not yet been done. 



2. (Ni\- x ,Pd x )2MnSn Heusler alloys 

The results of the pressure-dependence of T c for 
(Ni,Pd)2MnSn alloys are shown in Fig. [5] We have al- 
ready mentioned that the ratio by which T c increases 
with pressure for the end-point alloys is similar for both 
systems. This is also reflected in the concentration- 
dependence of T c , namely a similar increase of T c under 
pressure over the whole concentration range. The mono- 
tonic and essentially linear concentration-dependence of 
T c under pressure reflects a similar trend of exchange 
interactions found for the case of ambient pressure (see 
Refill). This similarity is due to the isoelectronic nature 
of Ni- and Pd-atoms coming from the same column, but 
different rows (3- and 4-rf, respectively) of the Periodic 
Table. 



3. (Ni\- x ,Cu x )iMnSn Heusler alloys 

Results of a similar study for (Ni,Cu)2MnSn alloy over 
a set of pressures realized by a linear reduction of the 
alloy lattice constant from the ambient value up to 3% 
are shown in Fig. [3] For comparison with experiments, 
results for the ambient pressure are slightly improved for 
the present LSDA+U model with respect to the LSDA 
results of our previous study.— However, both LSDA and 
LSDA+U reveal the same general trend. Specifically, 
we observe two concentration regions, the first one for 
x <0.4 in which T c is essentially constant with a shallow 



minimum and that for x >0.4 where T c increases mono- 
tonically with x . With increasing pressure the behavior 
in the first concentration region remains unchanged, but 
there is a gradual reduction of the slope in the other 
concentration region (the Cu-rich end). Such a develop- 
ment is clearly related to opposite trends in the pressure- 
dependence of T c in Ni 2 MnSn and Cu2MnSn, a reduction 
of T c with pressure for CU2M11S11 and an increase for 
Ni2MnSn. The crossover from the positive derivative of 
the pressure-dependence of T c to the negative one takes 
place around x=0.7. Thus, for the pressure correspond- 
ing to a 3% reduction of the lattice constant we predict 
only a weak dependence of T c on the alloy composition. 



C. Qualitative understanding of results 

In this subsection we will present a qualitative under- 
standing of the composition-dependence of T c at ambi- 
ent pressure as well as the effect of pressure on T c us- 
ing the idea of Anderson's superexchang o 27 ' 43 interaction 
and Stearns 4 ^ model of the indirect (RKKY-type) inter- 
action in ferromagnetic Heusler alloys. 

The Mn-Mn exchange interaction, as in most other 
cases, can be divided into three groups: direct, indirect 
and superexchangc. Because the Mn-atoms are not the 
nearest neighbors of themselves and have a large sep- 
aration (dMn-Mn > 4A), the direct interaction is not 
of any importance in Heusler alloys. The most impor- 
tant interaction turns out to be the RKKY-type indirect 
interaction^ Of somewhat minor importance is the su- 
perexchangc interaction, introduced in a series of papers 
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by Andcrson42 Superexchange, which is also an indi- 
rect interaction, is antiferromagnetic, while the RKKY- 
type interaction can be ferromagnetic or antiferromag- 
netic. Sasioglu et al£Z- recently used these two inter- 
actions to interpret the results for some Heusler alloys. 
Superexchange is related to possible transitions from oc- 
cupied states at the Fermi level Ep to unoccupied lev- 
els above, in the present case to the unoccupied minority 
Mn-bands. Because the superexchange interaction is neg- 
ative, the stronger the superexchange part, the smaller 
is the total exchange integral (and T c ). The strength of 
the superexchange interaction depends on the density of 
states (DOS) at Ep and the energy separation between 
Ep and the unoccupied Mn-bands. The larger the DOS 
at Ep and/or the smaller the above energy separation, 
the stronger is the superexchange. This lends a qualita- 
tive understanding of the increase of T c from Pd2MnSn 
to N^MnSn to Cu2MnSn, if the DOS mechanism dom- 
inates over the decreasing of energy separation in the 
above sequence of alloys (Fig.2 of Ref@). If we dope the 
Ni-sublatticcs with Cu-atoms in the N^MnSn host, the 
observed hockey-stick-like composition-dependence of T c 
can be understood as a competition of above two effects: 
for low Cu-content the decrease of T c with increasing 
Cu-concentration is due to reduction of the above en- 
ergy separation, causing an increase of the superexchange 
part of the interaction. For xcu > 0.35 the DOS-effect 
starts to dominate, with decreasing DOS at Ep and thus 
smaller superexchange and larger total exchange interac- 
tion. This leads to higher T c of Cu2MnSn as compared 
to Ni 2 MnSn. 

In (Ni,Pd) 2 MnSn alloy the DOS-effect is the only im- 
portant effect over the entire concentration range. Be- 
cause the Ni- and Pd-atoms have the same valency, the 
above mentioned energy separation does not change with 
concentration. There is a small decrease in the DOS 
at Ep in Pd 2 MnSn compared with Ni 2 MnSn (Fig.2 of 
Ref@), decreasing the superexchange and increasing T c . 

As mentioned above, although the superexchange 
mechanism seems to support the observed results, it can- 
not be the driving mechanism for the Heusler alloys under 
study. This is because of relatively large separation of the 
unoccupied minority Mn-peak above Ep (the Mn-DOS 
just above Ep is small). Of much-larger importance is the 
indirect interaction and its key feature in ferromagnetic 
Heusler alloys was captured and successfully described 
by Stearns^ in a series of papers in the late seventies. 
Stearns asserts that the magnetism of Heusler alloys (as 
well as bec Fe^ ) can be understood by dividing the d- 
electrons into localized d-electrons (di) associated with 
narrow c?-subbands and itinerant d-elcctrons (di) associ- 
ated with (usually) one broader d-subband. The inter- 
action between these localized and itinerant d-electrons 
is similar in nature to RKKY interaction, originally for- 
mulated for the interaction between localized moments 
in a free electron gas. The d-band gets narrower across 
(with increasing number of electrons) a given series and 
broader with increasing row number (e.g. from 3d to 



Ad). Thus Stearns used her model to describe correctly 
the change of T c from N12M11S11 to Ci^MnSn (an in- 
crease in di causing a rise of T c ) and from N^MnSn to 
Pd2MnSn (a decrease of T c due to decreasing d{). Uhl2£ 
used Stearns model to explain the hockey-stick appear- 
ance of the composition-dependence of T c with increasing 
x in (Nii-^jCu^MnSn Heusler alloys. With increas- 
ing Cu-conccntration, the lattice dilates. Initially, this 
dilation causes a reduction in the strength of the di-di 
interaction, associated with increased Mn-Mn distance. 
However, the addition of Cu also causes an increase in 
the number of di and thus in the di-di interaction. The 
resulting T c thus shows an initial decrease followed sub- 
sequently by an increase (for further details, see Uhl2& ). 
Thus both superexchange and indirect di-di exchange in- 
teraction give rise to the same trend for the composition- 
dependence of T c . 

The pressure effect on T c can be understood in the light 
of the above explanation for the composition-dependence 
at ambient pressure. Pressure reduces the lattice param- 
eter, an effect opposite to the dilation caused by Cu- 
addition. Pressure also broadens bands, reducing di and 
increasing di, again opposite to what happens with in- 
creasing Cu-content . Thus on the Ni-rich side increased 
pressure increases T c (opposite to the result with Cu- 
addition) , and beyond a certain critical Cu-concentration 
increasing pressure results in decreasing T c (opposite to 
the result with Cu- addition). 

Similar reasoning explains the composition- 
dependence of ambient pressure T c and the pressure- 
dependence of T c in (Nii-^Pd^MnSn alloys. Doping 
Ni-sublattice with Pd-atoms results in dilation of the 
lattice. In addition, d-electrons in Pd are more delocal- 
ized than in Ni. Both of these will lead to a decrease 
of T c with increasing Pd concentration. Pressure, via 
contraction of the lattice, reverses this effect. It appears 
that in this case the increasing strength of the interaction 
due to reduction in the Mn-Mn distance supersedes 
the effect of the pressure-induced derealization of the 
electrons. 

A couple of comments arc in order. Stearns model is 
applicable to only ferromagnetic Heusler alloys. In sev- 
eral scmi-Heusler alloys, or in general in situations where 
the superexchange part of the interaction becomes impor- 
tant, Stearns model may fail to predict the results cor- 
rectly. It would not be incorrect to point out that with- 
out an accurate first-principles study as presented here, 
it would be almost impossible to capture theoretically 
the details of the composition- or pressure-dependence of 
T c , although some qualitative understanding (especially, 
a posteriori) is possible. In this sense the importance of 
first-principles studies like the present one is irrefutable. 



D. Spin-disorder induced resistivity 

In magnetic alloys there are three contributions to 
the resistivity: a temperature-independent residual re- 
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r(K) 

FIG. 4: (Color online) The temperature-dependence of the re- 
sistivity of the (Nii-^Pd^MnSn alloys for ambient pressure 
(full symbols) and the pressure corresponding to the reduc- 
tion of the alloy lattice constant by 3% (open symbols): (a) 
x = 0, (b) x = 0.5, (c) x = l. 



sistivity due to the alloy disorder and other defects, and 
two temperature-dependent contributions, one due to 
electron-phonon and the other due to electron-magnon 
scatterings. In the present case, the residual resistivity is 
small because the states corresponding to Ni-, Cu-, and 
Pd-atoms are well below the alloy Fermi energy^ and 
thus influence states at the Fermi energy (relevant for 
residual resistivity) only weakly. The contribution due 
to phonons is well understood, and for Hcuslcr alloys 
around room temperature,— this contribution is known 
to be small. The dominant contribution is thus due to 
spin-disorder . 46 ' 47 The spin-disorder induced resistivity 
is controlled by the spin-spin correlation function which 
is small at and above T c , so that a reasonable repre- 
sentation of the situation is given by the DLM state, 
as explained in detail elsewhere* 6 - Kasuya's s-d interac- 
tion model^ 7 - shows that the temperature-dependence of 
spin-disorder resistivity is primarily quadratic and can 
be well-simulated by the £?T 2 -law. In the present case 
we determine the constant B from first principles using 
B = p(T c )/T^, where p(T c ) is identified with the resistiv- 
ity of the DLM state. The DLM state is described in the 
framework of the CPA as an equiconcentration alloy with 
two types of Mn atoms (A and B) with their moments 
pointing randomly up or down (collinear disorder). The 
corresponding resistivity can be properly calculated by 
the Kubo-Greenwood formula applied to this binary-Mn 
alloy. The results for Ni 2 MnSn, (Nio.5,Pdo. 5 )MnSn, and 
Pd2MnSn alloys are shown in Fig. [4] It should be noted 
that the residual resistivity of the (Nio.5,Pdo.s)MnSn al- 
loy in the ferromagnetic state is about two orders of mag- 
nitude smaller than that in the DLM stated 

The spin-disorder part of the resistivity is controlled 



by the constant B. For the cases shown in Fig. [4j T c 
increases under pressure. Due to band-broadening un- 
der pressure and accompanying derealization of the elec- 
trons, p{T) decreases in general for a given T less than 
the ambient pressure T c . However, p(T c ) (as well as p(T) 
for T > ambient pressure T c ) is higher under pressure 
because of the increase in T c . 

To our knowledge, there are no available experimental 
results for the systems studied in this paper, although 
the results of Austin and Mishra£ for related Pd2MnSb 
Heusler alloy do indicate, as in this work as well, an in- 
crease of the resistivity under pressure for temperatures 
above T c . 

Undoubtedly pressure has non-negligible influence on 
the part of the resistivity due to phonons. This, however, 
is beyond the scope of the present article. 



IV. CONCLUSIONS 

We have studied magnetic and transport properties of 
quaternary Heusler alloys (Ni,T)2MnSn (T=Cu, Pd) un- 
der pressure by means of the first-principles density func- 
tional method. In particular, we have investigated in de- 
tail the pressure-dependence of the Curie temperature. 
In agreement with experiments, we obtain an increase 
of the Curie temperature under pressure (reduction of 
the lattice constant) for Ni2MnSn and Pd2MnSn. On 
the other hand, Cu2MnSn alloy exhibits a reduction of 
the Curie temperature with applied pressure. The results 
can be qualitatively understood as an interplay of two ef- 
fects: an increase of the bare exchange integrals with the 
volume reduction and the decrease of magnetic moments 
with pressure due to band broadening. 

The concentration-dependence of the Curie tempera- 
ture under pressure is simple in (Ni,Pd)2MnSn Heusler 
alloys where the increase of the Curie temperature has 
an almost constant slope. On the other hand, we predict 
a dramatic change of the Curie temperature with pres- 
sure from the positive for Ni-rich (Ni,Pd)2MnSn alloys 
to a negative one in Cu-rich alloys, the crossover being 
around 70 % of Cu. We emphasize that such a com- 
plex behavior cannot be obtained without an accurate 
calculation of the exchange interaction and a proper sta- 
tistical treatment to compute T c (i.e., by going beyond 
the MFA). 

A simple explanation of the behavior of the calculated 
T c under pressure is given in the framework of Anderson's 
supcrexchange interaction and the Stearns model of the 
indirect exchange interaction between itinerant and lo- 
calized d-electrons. 

We have also investigated the pressure-dependence of 
the spin-disorder related part of the resistivity, which 
dominates at higher temperatures, the residual resistivity 
due to alloy disorder being small and important only for 
very low temperatures. We have found that the resistiv- 
ity at the Curie temperature increases with pressure, but 
this result is due to the increase of the Curie temperature 
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with pressure itself. For temperatures below the ambient 
pressure T c the resistivity decreases under pressure, as 
would be expected due to band broadening. 
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